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Objectives 


□ Support the Aeronautics Research Mission Directorate (ARMD) guidelines at NASA’s 
Dryden Flight Research Center. 

❖ Supported by Subsonic Fixed Wing (SFW) & Supersonics (SUP) projects under Fundamental 
Aeronautics (FA) Program 

□ Generating the basic object-oriented framework for a multi-disciplinary analysis and 
design optimization tool to be used in the preliminary design stage of a subsonic / 
transonic / supersonic / hypersonic aircraft. 

♦♦♦ Develop analysis modules for flutter optimization in transonic flight regime 

> No commercial MDO code exists in this speed regime 

♦♦♦ Difficulties in domain of analysis for each disciplines 

> Frequency-domain: Design, Classical control theories, Lifting surface theories, Linear, etc. 

> Time-domain: Analysis, Modem control theories, CFD, Nonlinear, etc. 

□ The framework will be set up to integrate analysis codes for multiple disciplines, instead 
of having one code perform the analysis for all the disciplines. These multiple analysis 
codes will be integrated using a front end code. 

□ Reduce uncertainties in the aeroservoelastic model to increase the safety of flight 

♦♦♦ Develop model update techniques based on design optimization to improve analysis/test 
correlation 
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Approach 


❖ 

♦>, 


❖ 


□ Front-End Code 

Based on FORTRAN “call system” command 
Optimizer 

> Genetic Algorithm: gradient free approaches 

> DOT : based on gradient values 

❖ Executive Pre- & Post-processors 

> Pre-processor 

■S Select optimizer, design variables, objective functions, 
and constraints 

> Post-processor 

■S Optimization histories & active constraints table 

□ Each Disciplines 

Pre-processor: Prepare input data for an analysis code 

> Read before optimization 
■S Analysis code name 

■S Temporary data file name created by other disciplines 
■S Code name for updating input data 

> Read during optimization 
■S Read updated values of design variables 
■S Update input data automatically based on new design 

Output 

> System variables (such as total weight, frequencies, drag, 
noise level, flutter speed, etc.) 

> Design sensitivity matrices if available 

> Temporary data file required by other disciplines 

Post-processor 

> Draw results using GUI 

> Create result data file 


Object-oriented MDAO Tool 




❖ 



Need to develop; 
Depends on 
Analysis code 


Analysis code; 
Commercial 
and/or in- 
house codes 
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Need to develop; 
Depends on 
Analysis code 

Chan-gi Pak-3 


Disciplines 


□ The following modules will be developed in the MDAO code. 

Flutter/Divergence analysis module 


❖ 


❖ 


❖ 


❖ 


❖ 


❖ 


Static structural analysis module 
Structural model update module 
Buckling analysis module 
Weight analysis module 
Modal analysis module 


❖ 


❖ 


❖ 


❖ 


❖ 


❖ 


❖ 


Structural 

Mode 

Interaction 

Test 


Trim analysis modul 
Sonic boom analysis 
Control model update module 
Aerodynamic load analysis module 
Aerodynamic model update module 
Gain/Phase margins analysis module 



Pre & Post 
processors 


Structural 

Model 

Update 



Ground 

Vibration 

Test 


Static Load 
/ Deflection 
Test 


'Deflection, 

: Stress, Strain, ] 
. & Buckling 


Frequency & 
Mode Shape 


Validated 
Steady & 
Unsteady 
Aero Model 





T 

Trim 

1 

L. 


J 


Validated 

Structural 

Model 


Flight 

Vibration 

Test 


Other performance analysis modules, such as cabin noise, mission analysis, landing and taxiing analysis, 
panel flutter analysis, hot structure dynamics, etc. 
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Analysis Codes for each Discipline 


Disciplines 

Analysis Codes 

FY 08 

FY 09 

FY10 

Stress/Strain Analysis 

NASTRAN sol 101 


In-house code 

Buckling Analysis 

NASTRAN sol 105 



Modal Analysis 

NASTRAN sol 103 


In-house code 

Lift, Drag, Trim Analysis 


FUN3D 

ZAERO, DLM, & 
KFM 

Subsonic/Supersonic Flutter Analysis 

ZAERO 

DLM 

KFM 

Subsonic/Supersonic Divergence Analysis 

In-house code 



Gain/Phase Margin Analysis 


In-house code 

ZAERO 

Transonic Flutter Analysis in Frequency-Domain 


In-house code 


Transonic Flutter Analysis in Time-Domain 


FUN3D 


Transonic Aeroservoelastic Analysis in Time-Domain 



Modify FUN3D 


Structural Dynamics Group 


Chan-gi Pak-5 


Validations 


□ 


□ 


□ 


Sample 1 : Weight minimization 
❖ Three bar truss 


Sample 2: Flutter/Divergence speed maximization 
❖ NASA 870 IKHANA Aircraft 


Sample 3: Minimize errors between test data and analytical results 

❖ Structural model update 

❖ Aerostructures test wing 1 
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Sample 1: Weight Minimization 


□ Three Bar Truss Problem 

♦♦♦ Check feasibility of our proposed approach 

♦♦♦ Test weight analysis module 

♦♦♦ Test static structural analysis module 

□ Objective Function: Minimize total weight 


□ Applied Load P 1 =P 2 = 20000 lb 


□ Constraints 

♦> Allowable stress 20000 psi (tension) 

♦> Allowable Stress -15000 psi (compression) 


□ Design Variables 

♦♦♦ Cross Sectional Areas X l = X 3 and X 2 
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Sample 1 : Results 


□ Closed form solution 

Assume: A, = X, 


2X I -2X 2 

V2(2X,-1) 


Objective f(X v X 2 ) = 2^2X l + X 2 

Function: 

2X, + 4lX 2 _ 

Active Constraint: ^=1 -> ^ = 


f(X l ) = 2y[2X l + 


2Xj - 2Xf 

V2(2A,-1) 


r V2(6X 1 2 -6X 1 +1 ) _ q 

1 (2X } -l) 2 


v _ 6 + Vl2 _ 

12 

X 2 = 0.408249 
/ = 2.638958 


0.788675 


□ Finite Element Model 

❖ Use MSC / Nastran SOL200 

♦> Use 3 CBAR (uniform beam) elements 
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Closed form 
solution 


MSC/ 

NASTRAN 


MDAO 
with DOT 


MDAO 
with GA 


Bar Area X 


.788675 


.77142 


.78798 


.80460 


Bar Area X 0 


.408249 


.45185 


.40999 


.36526 


Total Weight 


2.63896 


2.6338 


2.6388 


2.6419 


Number of 
iterations 


N/A 


50 generation 
& 

20 propulsion 


Three Bar Truss Optimization, Genetic, Stress 
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Sample 2: Flutter/Divergence Speed Maximization 


□ NASA Dryden Flight Research Center acquired a Predator B unmanned aircraft system for civilian 
missions. 


❖ IKHANA carries a ‘fire pod’ that will transmit images of remote areas of the western United States down 
from the aircraft to a ground station. 

❖ The fire pod is located under the wing near the left wing root, and can alter the flutter characteristics of the 
baseline aircraft. 


❖ Test modal analysis module. 

♦♦♦ Test flutter/divergence analysis module 


Baseline Configuration 




Structural Dynamics Group 


Chan-gi Pak-9 




Sample 2: Flutter/Divergence Speed Maximization (continued) 


□ 


□ 


□ 


□ 


Optimization Problem Statement 

♦> Objective Function: flutter and divergence speed 

❖ Constraint: None 

Design Variables 

❖ Chordwise location of the fire pod 

Structure Finite Element Model 

❖ MSC Nastran model (18854 nodes and 20979 elements) 

Unsteady Aerodynamic Model 

ZAERO model (2736 of elements) 


❖ 
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Sample 2: Challenges / Issues 


□ Modification of the fire pod location affects both the structural finite element model and 
the unsteady aerodynamic model. 

♦> New MSC/NASTRAN and ZAERO analyses must be executed for each optimization iteration. 
♦> Computing speed for generating AIC (from scratch 20 hours; even using existing 30 mins) 

♦♦♦ Genetic optimizer requires thousands of iterations 


□ Approximation Methods 

♦♦♦ Avoid computing a new AIC matrix for each design variable update. 

♦> AIC approximation based on matrix AJJ (General AIC). 

> ZAERO allows “Direct matrix input” for matrix QHH but not matrix AJJ. 

♦> AIC approximation based on matrix QHH (Modal AIC). 

> Cubic-spline each element in matrix QHH does not provide accurate estimates. (Mode switch etc.) 

❖ Flutter and divergence speed approximation based on pre-calculated values. 

> Interpolates flutter and divergence speeds from some pre-calculated flutter and divergence speeds for 
each design variable update. 
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Speed (Keas) Speed (Keas) 


Sample 2: Results 


Flutter & Divergence Speed 



Divergence Speed 
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Ikhana with fire pod optimization using 
design optimization tools. 


Design variable 

Chordwise location from baseline, in. 

Objective value 
Critical speed, KEAS 

Initial value 

Final value 

-20.0 

-19.80 

463.95 

0.0 

0.0001 

464.20 

10.0 

14.31 

465.53 

20.0 

14.27 

465.53 


Summaries of critical speeds before and after optimization 
using genetic algorithm. 



Design variable 
Chordwise location 
from baseline, in. 

Objective value 
Critical speed, KEAS 

Baseline 

0.00 

464.20 

Optimized 

14.37 

465.53 
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Sample 3 : Structural Model Update 


□ Everyone believes the test data except for the experimentalist, and no one believes the 
finite element model except for the analyst. 

♦♦♦ Some of the discrepancies come from analytical Finite Element modeling uncertainties, noise 
in the test results, and/or inadequate sensor and actuator locations. 

□ MIL-STD- 1 540C Section 6.2.10 

❖ Test Requirements for Launch, Upper- Stage, & Space Vehicles 

♦♦♦ Less than 3% and 10% frequency errors for the primary and secondary modes, respectively 

♦> Less than W% off-diagonal terms in orthonormalized mass matrix 

□ AFFTC-TIH-90-00 1 (Structures Flight Test Handbook) 

♦♦♦ If measured mode shapes are going to be associated with a finite element model of the 
structure, it will probably need_ to be adjusted to match the lumped mass modeling of the 

analysis. 

♦> Based on the measured mode shape matrix [F] and the analytical mass matrix [M] , the 
following operation is performed. 

0> T M<5 

♦♦♦ The results is near diagonalization of the resulting matrix with values close to 1 on the 
diagonal and values close to zero in the off-diagonal terms. Experimental reality dictates that 
the data will not produce exact unity or null values, so W percent of these targets are accepted 
as good orthogonality and the data can be confidently correlated with the finite element model. 
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Sample 3 : Structural Model Update (continued) 


□ Flutter Analysis 


❖ 


❖ 


Uncertainties in the structural dynamic model are eliminated 
through the use of “model tuning technique” 


Flutter Analysis Procedure 
(ax NASA DFRC 


❖ 


Based on analytical modes 

□ Update Finite Element Structural Model using Test Data 

Recall MIL-STD- 1 540C Section 6.2.10 

> Less than 3% frequency error: primary modes 

> Less than 10% frequency error: secondary modes 

> Less than 10% off-diagonal terms in mass matrix 

Update Mass 


Measure Weight, C.G., & 
Moment of inertia 


❖ 


GVT 


O 


G> 


(Or 


❖ 


> Minimize errors in total weight, C.G. location, and mass moment 
of inertia 

> Minimize off-diagonal terms in orthogonal mass matrix 

Update Stiffness 

> Minimize errors in frequencies 

> Minimize errors in mode shapes and/or minimize off diagonal 
terms in orthogonal stiffness matrix 


Wind 



C 


FE Model M& K 


Start 


J> 


Optimization Step 1 : 
Update Mass Properties 

W, X CG , Y cg , Z cg , Ixx, 
Iyy, Izz, Ixy, Iyz, & Izx 

I 


Optimization Step 2: 
Improve Orthogonality 
of Mass Matrix M 


Optimization Step 3: 
Update Frequencies 
& Mode Shapes 


C 


End 


J> 


Aero Model A 

~~y~ _ 


Flutter Analysis 
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Sample 3: Mathematical Background 


□ Optimization Problem Statement 
Minimize J 


❖ 


❖ 


Such that 


h — £ k 


k = 1 ... 1 3 & k^i 


□ Step 1 : Improve Rigid Body Mass Properties 
Errors in Total Mass 
Errors in CG Locations 
Errors in Mass Moment of Inertias 


❖ 


❖ 


❖ 


Mass Properties 

Objective Functions & Constraints 

Total Mass 

/, = (W-W g ) 2 /W g 2 

CG Locations 

J 2 = (X-X g ) 2 /X g 2 

j 3 = (Y-Y 0 )W 0 2 

J 4 = (Z-Z 0 ) 2 /Z 0 2 

Mass Moment of Inertias 

j 5 = (Ixx-Ixxg) 2 /Ixxo 2 

J 6 = (Iyy-Iyyg) 2 /Iyyg 2 

— (^ZZ - ^ZZg) 2 ^ZZG 2 

*^8 _ (^XY"-^XYg) 2 ^XYG 2 

J9 — (^YZ - ^YZg) 2 ^YZG 2 

J 10 ~ (^ZX"-^ZXg) 2 ^ZXG 2 
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Sample 3: Mathematical Background (Continued) 


□ Step 2: Improve Mass Matrix 

♦♦♦ Off-diagonal terms of Orthonormalized Mass Matrix: M = <D G T T T M TO G 
Guyan reduction 



I 


\M 

M 1 


[K 

K 1 

II 

ii 

— K^K 

^ Vs/;/ 

M = 

mm 

ms 

K = 

mm 

ms 

yj 


M sm 

M ss 


TS 

V sm 

1 


Improved reduction system 


T = T = 

1 1 IRS 


I 

-K;X m +(K;:M ss 


g 1 K g 


M g = T g MT G 


K g =T t g KT g 


j„= Em- 



FEM 



GVT 
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Sample 3: Mathematical Background (Continued) 


□ Step 3: Frequencies and Mode Shapes 

Errors in Frequencies 


❖ 


n 


J, Z = X 




i = 1 


i i 


V 


Q 


y 


❖ 


Option 1 : Off-diagonal terms of Orthonormalized Stiffness Matrix: K= O g t T T K TO 


G 


n 


j„= X K 

i=\,j=\,i^j 


Option 2: Errors in Mode Shapes 


m 


j n =y,(®,-' k) 


i — 1 


n: number of modes m: number of sensors 
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Sample 3: Aerostructures Test Wing 1 


Number of DOFs in FE Model 

1311 

Number of Accelerometers for GVT 

35 


AEROSTRUCTURES TEST 

WING 


Drydsn 


April 24, 2001 


cc 




NASTRAN Structure Model 
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Sample 3: Results 


Frequency Comparisons 



GVT 

Before Optimization (ATA’s final) 

After Optimization without ./, j 

After Optimization with J x , 

Frequency 

(Hz) 

Frequency/Error (Hz/%) 

MAC 

Value 

Frequency/Error (Hz/%) 

MAC 

Value 

Frequency/Error (Hz/%) 

MAC 

Value 

Guyan 

Full order 

Guyan 

Full order 

Guyan 

Full order 

Mode 1 

13.76 

13.35/-3.0 

13.35/-3.0 

99 

13.75/-0.1 

13.75/-0.1 

98 

13.41/-2.5 

13.41/-2.5 

95 

Mode 2 

20.76 

22.82/9.9 

22.82/9.9 

99 

20.76/0.0 

20.76/0.0 

99 

21.01/1.2 

21.01/1.2 

97 

Mode 3 

77.83 

79.06/1.6 

78.77/1.2 

95 

77.82/-0.0 

77.84/0.0 

95 

77.87/0.1 

77.50/-0.4 

95 


Mode 1 



L 

First Vertical Bending 



Mode 2 


Second Vertical Bending 


1.50-01 l 


-1.00-03 



Mode 3 



First Torsion 
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Sample 3: Results (continued) 


Mass Properties 



Measured 

Before Optimization 

After Optimization without J x x 

After Optimization with J x x 

Weight 

2.66 lb 

2.77 lb (error 4.1%) 

2.67 lb (error 0.4%) 

2.70 lb (error 1.5%) 

x CG 

N/A 

12.94 inch 

12.88 inch 

12.72 inch 

y cg 

N/A 

9.16 inch 

8.80 inch 

8.91 inch 

7 

^CG 

N/A 

0.0 inch 

0.0 inch 

0.0 inch 

^xx 

N/A 

161.22 lb-inch 2 

152.06 lb-inch 2 

154.78 lb-inch 2 

^YX 

N/A 

95.27 lb-inch 2 

93.75 lb-inch 2 

89.45 lb-inch 2 

I yy 

N/A 

1 13.08 lb-inch 2 

112.83 lb-inch 2 

102.57 lb-inch 2 

^zx 

N/A 

0.011 lb-inch 2 

0.010 lb-inch 2 

0.010 lb-inch 2 

^ZY 

N/A 

-0.028 lb-inch 2 

-0.035 lb-inch 2 

-0.033 lb-inch 2 

^ZZ 

N/A 

268.2 lb-inch 2 

258.79 lb-inch 2 

25 1 .26 lb-inch 2 

Orthonormalized 
Mass Matrix 

1 

8.9 % 

17.7 % 

1 

15.7% 

14.8 % 

1 

3.0% 

11.9% 

.089 

1 

9.3 % 

.157 

1 

10.9 % 

.030 

1 

2.8 % 

.177 

.093 

1 

.148 

.109 

1 

.119 

.028 

1 
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Conclusions 


□ An object-oriented MDAO tool has been developed at the NASA Dry den Flight Research 
Center to automate the design and analysis process and leverage existing commercial as 
well as in-house codes to enable true multidisciplinary optimization in the preliminary 
design stage of subsonic, transonic, supersonic and hypersonic aircraft. 

♦> Once the structural analysis discipline is finalized and integrated completely into the MDAO 
process, other disciplines such as aerodynamics and flight controls will be integrated as well. 


□ Simple and efficient model tuning capabilities based on optimization problem are 
successfully integrated with the MDAO tool. 


♦> More synchronized all phases of experimental testing (ground and flight), analytical model 
updating, high-fidelity simulations for model validation, and integrated design may result in 
reduction of uncertainties in the aeroservoelastic model and increase the flight safety. 
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Questions ? 



Pre & Post 
processors 



Structural 

Model 


Ground 

Vibration 

Test 


Gain/Phase 

Margin 


Structural 

Model 

Update 


r i 

Static Load 


/ Deflection 

L J 

Test 


^Aerodynamic^ 
Model 
^S^Update^ 

4 A k 


Flutter/ 

Divergence 


Validated 
Steady & 
Unsteady 
Aero Model 


Lift & Drag 


'Deflection, 
Stress, Strain, ] 
. & Buckling 


Frequency & 
Mode Shape 


'Aerodynamic 
Model 


Validated 

Structural 

Model 


Flight 

Vibration 

Test 
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